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TRAVELING-WAVE CHARGED PARTICLE 
GENERATION 



LANCE HEROLD 



ABSTRACT 



Submitted, to the Department of Naval Architecture and 
Marine Engineering on May 19, 1967 in partial fulfillment 
of the requirements for the Master of Science Degree in 
Electrical Engineering and the Professional Degree, Naval 
Engineer. 

A sinusoidally charged aerosol is injected into a fluid 
flow and interacts synchronously with a traveling wave of 
voltage. The voltage traveling wave can either be that 
induced on. passively loaded channel electrodes, or externally 
supplied and adjusted to either lead or lag the charge wave. 

In both the passive loading and lagging impressed-voltage 
cases, fluid dynamic to electrical energy conversion is 
achieved, while in the case of voltage wave leading the 
charge wave , electrical pumping of the fluid flow occurs. 
This study concentrates on an induction method of charging 
a distilled water aerosol formed in a set of spray atomization 
nozzles using compressed air. 

Historical background of direct-current electro-fluid 
dynamic energy conversion is outlined, along with a 
theoretical comparison with what is considered a hitherto 
untried synchronous approach to the problem. 
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Experiments were conducted to determine a feasible 
method of inducing alternating charge on aerosol droplets; 
small enough to produce useful concentrations, and large 
enough for substantial viscous coupling with the transport 
medium. The resulting charged aerosol generator is described 
and analyzed. 



Thesis Supervisor: James R. Melcher 

Title: Associate Frofessor of Electrical Engineering 
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I . INTRODUCTION 



A • Back gr ound 

Mechanical to electrical energy conversion systems have 
traditionally employed the interaction of charge in motion 
and magnetic fields. Since magnetic field energy densities 
are characteristically orders of magnitude greater than those 
of electric fields, relatively little attention has been 
directed toward electric field energy conversion machines. 

A characteristic of Electro-Fluid Dynamic systems in 
general, has been their high voltages (tens of kilovolts) and 
small currents (microamp to milliamp regime)- Electrostatic 
effects are often observed in corona discharge from long wires 
during snow storms, lightning, and other natural phenomena which 
seem to have inspired some of the early work in electric field 
energy conversion. • The earliest converters were Van de Graaff - 
type systems, in that they employed frictional or corona 
electrification and mechanical transport of the charged particles 
to a region of high potential. One of the first EFD generators 
was proposed by R.E. Vollrath^^ in 1932. His was a high 
voltage generator using blown powder, contact electrified, and 
carried to a metal sphere at high electric potential. This 
model v/as more in line v/ith present work, since it involved 
gas dynamic propulsion of the charges. 

The literature contains few, if any, significant develop- 
ments in the field from 1932 until the early 1950 s. In 1953 > 



2 



(n) 

A. M. Marks v J obtained, a patent for a thermo-electric power 
converter using a charged aerosol. Others had begun taking 
at least a limited second look at EFD machines as a possible 
method of producing electrical power without mechanical 
complexity, and hofefully, at higher power to weight ratios 
than are obtainable with magnetohydrodynamic systems. 

Within the last decade there has been a burst of activity 
in the study of electro-fluid dynamic energy conversion processes. 
The impetus for such research has come largely from the need 
for simple, light weight and reliable generators to be carried 
in space vehicles. An EPD-type system would uniquely meet 
these particular requirements if sufficient power levels and 
electrical efficiencies could be achieved. To this end, a 
number of feasibility studies have been conducted and various 
types of experimental systems have been constructed. It 
should be noted that all published, work, to date, seems to 
have been limited to uni-polar (d-c) EFD systems. 

B. Feasibility Considerations 

1. Mobility 

The most basic limitation on EFD energy converters 
is their extremely low efficiency in terms^ of electrical energy 
output for fluid dynamic energy input. A major reason for this 
is the limited nature of the interaction mechanism between 
charge carriers and the transporting medium. This has 
historically been expressed as mobility (K) , defined as the 
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ratio of drift velocity to the electric field. 

u 

K - (1.1) 

E 

Since the electric force on a charge carrying particle is 
the charge multiplied by the electric field, 

F = qE (1.2) 

the mobility can be expressed as: 

uq 

K = (1.3) 

F 

The nature of mobility for various fluid working media and 
particle sizes is described in several of the references and 
summarized in Appendix A, The facts which emerge from mobility 
consideration are that energy conversion efficiency increases 
as mobility decreases, and that the best coupling forces are 
achieved when the particles are large enough to obtain viscous 
interaction. The latter has been found to require particle 
sizes greater than about 1 micron diameter^’ . 

Using Stokes' lav;, the force balance equation for a 
particle becomes: 



F = qE = 6xuap. 



(1.4) 
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Combining equations (1.3) and (1.4): 

q 

K « (1.5) 

6 tx y. a 

The fact that mobility varies inversely with viscosity 

led to some experimentation with non-conducting liquids. In 

(3) 

19&1» Stuetzer w/ presented the results of theoretical and 
experimental work v/ith a small high voltage generator using 
transformer oil, kerosene, and air as the working media, and 
ions as charge carriers. 

The use of liquid transport media has two disadvantages 
which seem to have limited its consideration, especially for 
space vehicle generators or propulsion systems. These are the 
extra added weight and the increased wall friction losses. 

There are, in general, two ways of minimizing the effects 
of mobility in a gaseous medium. 

(1) Increasing operating pressure and flow velocity 
to decrease mobility (see Appendix A), and make 
drift velocities negligible v/ith respect to the 

• flow velocity. 

(2) Increase droplet size (at constant charge per 
droplet) to reduce mobility in accordance with 
Stokes* lav/. 

In the case where ions are the charge carriers, method (1) 
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is usually applied ^ A combination of the two 

methods is contained in the device described by Marks, Barreto, 
and Chu in Reference 7 * As might be expected, however, 
factors which decrease mobility often have an adverse effect 
on other parameters of the system. 

2. The Ballistic Approach 

Mobility can be reduced or removed entirely as a 
factor in system performance by converting only the kinetic 
energy of the charge carriers into electrical energy. This 
is accomplished by injecting colloidal or ionic charged particles, 
at high velocity, into a very low pressure converter where 
viscous effects would be minimized. 

Hasinger^ 6 ^ gives a detailed theoretical analysis of 
this " Ballistic^ type system, where he represents the energy 
conversion by equating the electrical power output to the 
change in kinetic energy between inlet and exit electrodes. 

4 V = i m (vf - vf) (1.6) 

where : V = electrical potential between inlet and 

outlet 

He also points out that an inherent loss comes from the 
necessity of allowing sufficient exit velocity to overcome 
space charge effects at the electrode, and variations in 
particle sizes. Hasinger concluded that 25 per cent of the 
kinetic energy should remain at the collector to meet these 
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considerations. Accordingly, v 2 v/ill be half of 

By considering space charge density in the conversion 
section to produce an electric field exactly offsetting that 
of the collector, Hasinger arrives at a simple expression for 
the upper limit on the electric force concentration in the 
ballistic generator. 

e f 2 

Force/unit area = - 1 (1.7) 

2 

where: E^ is the field at the entrance electrode, 

acting to decelerate the charged particles 



With the 25 percent loss, this leads to an expression for 
power concentration 

• Power/ unit area = E7 U / T o\ 

g 1 o vi.o; 

for v/hich Hasinger concludes that concentration of 300 
2 

watts/cm should be feasible. 

A type of ballistic EFD generator is described in 
Reference 6 and further discussed in Reference 8 . Disadvantages 
characteristic of ballistic systems are the requirement for 
focusing the particle beam, and difficulties inherent in the 
basic concept for recirculation or successive staging in a 
closed, low pressure system. 



3 . Vis cous-Coupl e d System Performance 
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The basic theoretical performance capability of an 
EFD system is derived in Appendix A. 3 following the simplified 
model used by Lawson The resulting expression for the 

maximum power per unit area which can ideally be extracted 
from the flow, in the form of pressure drop times velocity, 
is 

e ° a 

a H 0 =- e2-U 0 (1.9) 

This result is not surprising, but indicates the nature of the 
problem which presents itself when one contemplates the 
extraction of significant electrical energy from fluid dynamic 
head. As Lawson observes, when the parameter values corresponding 

to atmospheric pressure are inserted into equation (1.9), 

-4 

the maximum pressure drop is about 4 x 10 atmospheres. 

By increasing the pressure to 30 atmospheres the maximum &p 
increases to about one third of the standard atmosphere. 

This type of analysis seems to indicate, therefore, that in 
order for such a system to attain practicability, some kind 
of multi-staging in a closed system would be required. 

The concept of multi-staging imposes additional limitations 
on system performance. Since frictional power losses are 
proportional to (U )^, each collector grid and attractor 
electrode added to the channel further restricts the allowable 
flow velocity. 
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Recognizing this, Lawson defines stage efficiency as 



where : L 



'st 



L f + L d 



(I. 10) 



is electrical power density (see equation 
(A. 28)) 

L^. is fluid power input (1.9) 
is fluid dynamic power loss 



U: 



L d = “ Po 



(I. 11) 



where: a is a constant coefficient appropriate for 

the particular electrode and channel wall 
configuration. 

The expression for the maximum attainable stage efficiency 
is then, 



r 4 < u o - KE b> 



l st 






u? 



U + a 



m 



or 



K 



1 - 



Eb 



U. 



'st 



l + a P 






e o 4 



( 1 . 12 ) 
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Using air at 30 atmospheres pressure as the working medium, 
with flow speed of 50 m,/sec., and colloids as charge carriers 
for negligible mobility, Lawson arrives at a representative 
stage efficiency of 80 percent with a power concentration 

p 

of 200 watts/cm . 

It must be emphasized that this is a theoretical limit, 
does not consider exit velocity as a loss, and that present 
state of the art is a long way from physical realization of 
this performance. 

0 . P ROPOSED SYSTE?1 

Operation of an EFD energy conversion system for "alternating 

current " charged particles has been virtually ignored. Von 

( 12 ) 

Ohain and Wattendorf recommend such research^ ' , but 
apparently only within the context of the system described in 
the previous section. That is, the same transport and 
collection mechanisms would be employed. 

The synchronous EFD system investigated and described 
in this report, and in Reference 19, is based on an induction 
concept which is believed to be unique in this application, 
and possessed of distinct advantages over existing systems. 

1 . Descriptio n 

Sinusoidally charged droplets are injected into a 
fluid dynamic flow to create a traveling wave of charge confined 
in a channel. Electrodes are suitably placed along the 
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channel walls to interact synchronously with the charge 
wave; either to convert fluid dynamic to electrical energy, 
or vice-versa, in which case the system acts as a pump. 

There are basically three " modes " of operation for such 
a system. 

1. Channel electrodes are passively loaded (R,L,C) 

2. . A traveling wave of voltage is impressed on the 
electrodes which leads* the wave of charge 
3. The impressed voltage wave "lags” the wave of 

aV n — A 

V • t 

The basic theory of the energy conversion mechanism is 
derived in Appendix A and discussed in greater detail in 
Reference 19. The form of the solution (equation A. 19, page 4-8 ) 
suggests the analogy betv/een this synchronous EFD system and 
the magnetic field/torque angle relationships which obtain 
in conventional rotating machinery systems. 

In this regard, modes (1) and (2) correspond to a 
synchronous motor, while mode (3) corresponds to a synchronous 
generator. The channel electrodes are analogous to the stator 
and the fluid flow corresponds to the rotor. 

2 . Advan ta ges of Synchronous System 

The fundamental reason for interest in alternating 
current EFD power generation is the possibility it offers for 
obtaining a-c power in a system requiring neither moving 
parts nor electronic conversion of d-c. In addition to the 
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previously mentioned advantages of the direct EFD systems, 
the synchronous generator concept inherently resolves several 
of the major difficulties. 

Because the charges need not be collected in order to 
produce electrical power, the channel can be kept free of 
obstructions downstream from the charge wave source, and the 
system frictional losses are produced by wall friction alone. 
Multi-staging is a built-in feature of the concept, since 
each wave length of electrodes is equivalent to a stage. 

The principal restriction on flow velocity comes from the 
relationship between velocity, frequency, and wave length in 
a traveling v/ave. Thus, at a given frequency ,• tne length of 
channel per stage is proportional to U 0 « 

The flexibility of this system, in its capability for 
pumping non-conducting fluids as well as for electrical power 
generation, may eventually prove to be a significant advantage 
Another restriction on synchronous system performance 
comes from lateral stress on the charge flow due to the 
wall electrodes. This is in addition to the space charge 
forces existent in the d-c converter, and may prove to be a 
restriction on the number of wave lengths which can be accommo 
dated between renewals of the charge v/ave. 

The intention of this thesis research was to experi- 
mentally test the conversion mechanism theoretically predicted 
and to investigate an induction method of generating the 
charged droplet wave. 



II. EXPERIMENTAL PROCEDURE 



A . Description of Initial Exner i m e nts 

Initial experimentation was conducted using droplets 
issuing from capillary tubes placed directly in the air stream 
from the axial fan. At this stage the objective was to use 
the maximum droplet size which did not result in wetting of 
the channel walls and consequent shorting of the output 
electrodes. 

Various means of introducing the water droplets and 
various inducer electrode configurations were investigated. 
These included from one to five - 10 mil inside diameter 
capillary tubes mounted in a 1/8 inch O.D. tube which v/as 
inserted either from the side of the channel or from the narrow 
ends. Because of turbulence created in the channel, the latter 
method proved more satisfactory from the standpoint of channel 
wetting. For these experiments the excitation was provided 
by a 2 to 5 kilovolt, 60-cycle signal applied to a pair of 
wires in mid-channel, running vertically one centimeter apart, 
and about .6 cm. from the ends of the tubes. The experimental 
apparatus is shown in Figure 1. 

The droplets were drawn out of the capillary tubes by 
a combination of hydrostatic head, eductor action of the 
air stream, and electrostatic stress as described in Appendix 
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B. Of these, the hydrostatic head was by far the major 
factor, and variations in output current were observed to 
conform to the rate of droplet formation due to variations 
in v/ater pressure. It was also readily apparent that the 
maximum charge per droplet was at least an order of magnitude 
less than Rayleigh's limit and imposed no restriction on 
the results. 




Figure 1 . CAPILLARY TUBE AND INDUCER APPARATUS 

For LARGE droplets 
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The droplets were collected at the end of the channel 
by a wire mesh suspended from insulators within a shielded 
box and discharging through an impedance of ,65 megohms. 

Crude as the instrumentation proved to be, the quality of the 
measured signal could be observedly improved by increasing 
the number of droplets, since the rate of emission of droplets 
must clearly be greater than 60 per second if a relatively 
smooth 60 cycle current is to be produced. 

With droplets as large as those emerging from the smallest 
capillary tubes, any contact with the channel walls is 
unacceptable. This aspect of the problem proved go be such 
a dominant constraint that it dictated all subsequent designs. 
As an illustration of the effect of turbulence, a 28 mil 
0,D, capillary tube was mounted in the end of a one eighth inch 
feed tube which was then inserted into the 60 feet per second 
air flow, just downstream of the flow straightener section 
in the 2 1/4 inch by 6 inch channel. The feed tube entered 
from the middle of the 6 inch, '*near"side of the channel. 

With the mouth of the capillary tube positioned at mid- 
channel there was no wetting of either side in the 36 inch 
downstream electrode section. The turbulence at the near wall 
prevented the droplets from being positioned more than about 
1/8 inch off-center in that direction without wetting in the 
channel. The absence of additional turbulence at the far 
wall permitted movement of the droplet emitter to within 
about 3/8 inch of that side before the onset of wetting 
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down channel . 

Because of this limitation on the distribution of large 
droplets no further work with them was attempted in the course 
of this particular project. This is not to say that such 
a study would necessarily be impractical, however, and 
further commentary on this subject is made in the recommenda- 
tions of Chapter VI. 

B , Small Dro plet Exp eriments 

Considerations 

As described in section C of Chapter I, mobilities 
deriving from the application of Stokes' law are achieved 
as long as droplet diameters are on the order of one micron 
or greater. Since a complete investigation of system per- 
formance using the entire range of droplet sizes v/as not 
feasible for this study, it v/as decided to concentrate on 
the production of charged “fog /( droplets ^ approximately 1 to 
20 microns in diameter. The objective was to fill the channel 
with droplets large enough to obtain satisfactory viscous 
coupling with the air stream, but small enough to permit 
evaporation to offset the deposition of droplets on the 
plexiglass wall surfaces. 

Another advantage of using fog droplets comes from the 
signal smoothing effect of a high rate of charge carrier 
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formation. An obvious disadvantage is the greater difficulty 
of inducing alternating charge on aerosol sprayed droplets. 
Learning the extent of this problem and the possibilities of 
solving it became the next experimental objective. 

Spray Nozzl e Experiments 

A commercial type spray nozzle was used for feasibility 
investigation. The particular configuration for air and water 
mixing is shown in Figure 2. It was assumed that this method 
of fog generation would provide the best type of water jet 
break up from the standpoint of access -i >>-5 litv to the inducing 
f 5 eld , 

Nozzles which mix air and water before emission could be 
expected to be less desirable because of the poor conductivity 
of the fog, even at the initial concentration. Another method 
of aerosol generation consists of intersecting air and water 
jets, but it was not tried for this narrow channel because of 
the better directional control afforded by a nozzle. 

The other advantageous feature of the concentric air jet 
is the shield it provides around the water stream. This 
enables the inducer electrode to be positioned close to the 
region of the fog formation with minimum possibility of an 
electrical short circuit between electrode and water. This 
was tested by mounting the commercial nozzle in a piece of 
plexiglass with a circular inducing electrode of 3/16 inch 
diameter centered on the nozzle outside of a 1/16 inch 
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plexiglass separation. 




Experiments did show, that significant charge can be induced 
on the fog droplets in this way. 

The results of the experiments with the commercial nozzle 
are plotted in the next chapter in Figure J. The current 
produced by this configuration was about 2 orders of magnitude 
greater than that produced by the capillary tube experiments. 
It should be noted that the measurements of Figure 3 were 
made with the nozzle outside of the channel and are the total 
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current produced by one nozzle. A rough drawing of this 
nozzle and electrode is presented in Figure 10. 

Nozzle Design for the C h arge Wave Gener at or 

Having established the feasibility of using aerosol 
spray technique, the next step was to achieve a design v/hich 
could fill the channel v/ith droplets with minimum interference 
with the air stream and v/etting of the channel walls. Catalog 
information for the commercial nozzle describes the qualitative 
effect of variations in the air pressure with respect to 
flow rate and average droplet size. The effect of increasing 
air flow, at constant water siphon height, is to retain 
or slightly decrease the droplet diameter while increasing 
the water flow rate. 

Several different nozzle designs were constructed using 
capillary tubes of varying diameter together with various 
sized air annuli. The process leading to the result shown 
in Figure l6 was largely trial and error, first using steel 
and plastic tubing, and finally plexiglass and stainless steel 
capillary tubes. 

Two different sets of plexiglass and capillary tube 
nozzles finally resulted. The first consisted of five 
vertically aligned nozzles in each of two 3/16 inch thickness 
plexiglass sections. These used 10 mil inside diameter 
capillary tubes and were designed to produce a fine spray 
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on the order of a few microns diameter. The second was similar, 
"but contained 10 nozzles in each bank and used 15 mil inside 
diameter tubes to produce a great quantity of larger droplets, 
on the order of 20 microns diameter. 

Another feature which distinguished the two designs was 
the inducing electrode configuration. The fine droplet design 
used a flush, painted electrode surrounding the air ring , 
while the other consisted of 72 mil O.D. tubing sections, 
of about 2 nun length, concentric to the nozzle aperture. 

A • tailed description of the nozzle and electrode designs, 
with diagrams of the apparatus, is conoaxiicd in Appendix B.2. 
i'notographs of the generator apparatus, 20-nozzle unit, and 
generated aerosol, are included in the next chapter. 



III. RE SULTS 



A . General Description 

The charged aerosol generator units which finally evolved 
in the course of experimentation were successful in producing 
a significant charge wave of droplets of the desired size. 

In addition, the input impedance is infinite as frequency 
goes to zero, consisting of n pure" capacitance, as v/as desired. 

In order to evaluate the generator design it was necessary 
to make a number of measurements of the system parameters. 
Because the nozzles and electrode S VV \J V Zs xj-ckaI it was 

difficult to get completely consistent measurements from day 
to day. Non-unif ormities such as capillary tubes being 
slightly off-center in the air jets and electrodes, variations 
in the pressure from the air supply, and inherent difficulties 
in measuring small currents in an aerosol spray over a range 
of flow conditions, were some of the contributing factors. 

The data summarized in Appendix D is, for the most part, 
the result of 4 or 5 reasonably consistent test series. Since 
the objective of the project was primarily exploratory in 
nature; to investigate feasibility of the general concept, 
sophisticated instrumentation v/as not employed, or considered 
a productive goal at this stage of system development. 

Photographs of the apparatus and generated aerosol are 
displayed in Figure 5» 
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B ♦ Graphi c al Summary of Res u lts 
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SYNCHRONOUS GENERATOR APPARATUS 





CHARGED AEROSOL Generated 
by 10-nozzle unit 



Figure 5 



20-NOZZLE CHARGED 
AEROSOL UNIT 
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C . Elec tr ical Equivalent Circuit 

Ideally, the input terminal characteristic of the charge 
wave generator should be a "pure” capacitance, all of which 
is attributable to the fluid between inducer electrodes 
and aerosol droplets in the process of formation. In this 
ideal case, input current from the voltage source would equal 
the output current from the generator. 

Infact, there is also capacitance between electrodes and 
capillary tubes, and between leads. That fraction of input 
current attributable to this portion of the total admittance, 
and therefore unproductive of output """"wt. must be considered 
a system loss. 

Minization of the undesirable capacitance will be an 
important consideration in realization of an over-all system 
capable of self excitation. There are, therefore, essentially 
two electrical characteristics of the generator indicative of 
its performance. These are: 

(1) The ratio of current output to total input 
current s q- 

(2) The ratio of current output to voltage input 

= P 

In the case of the lOnozzle aerosol unit, at a frequency 
of 100 cps, 40 psi air, and input voltage of 65 v. r. m, s.' , 
input current v/as 1.16 panes with output current of .75 uanp. 
This results in a value of t\ i = 64.5 percent. ? for this unit 



\ 



PHASE ANGLE 
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Figure 6. BODE PLOT OF FREQUENCY RESPONSE 
for 10-nozzle unit 
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is equal to the slope of the I vs Vh curve in Figure 3» 

The result is J3 = ,ol2 paiap/volt, where current is measured 
peak-to-peak and voltage is root-mean square. 

Corresponding values for the 20 nozzle unit, using input 
voltage of 120 v-rms and 20 psi air, were: = 1,65/2.5 = 66 

percent, and 3 = .014- pamp/volt. 

Equivalent Capacitance 

Details of the measurements and calculations leading to 
an estimate of the input capacitance are contained in Appendix 

C. 

»« I* 

Figure 6 shows a Bode plot of the frequency response 
data, summarized in Appendix D, Using the order of magnitude 
capacitance predicted by the theoretical model of Appendix C, 
a value of series resistance was selected to place the break- 
point within the range of available frequencies. For an 

equivalent length of about 1/2 cm, equation (C,9) predicts 

-12 

a capacitance of about 4 x 10 farad. Accordingly, a series 

resistance of 48.5 meg-ohms was used in the measurement 

apparatus of Figure 18. The equivalent capacitance indicated 

-12 

by Figure 6 is C = 7*5 x 10 farad. 

Direct calculation of capacitance using equation (C.5) 

for a frequency of 125 cps and series resistance of 10.2 

-12 

meg-ohms, yielded a value, C ^ 6.65 x 10 . 

The value which will be taken as the average of these 

-12 

results, for the 10 nozzle unit, is C = 7 x 10 farad. 
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The equivalent model for the complete charged droplet 
generator will be the simple voltage-controlled current 
source shown in Figure 7* 




Figu re . . 7 . CIRCUIT MODEL FOR CHARGED DROPLET GENERATOR 



where C = 7 x 10 -1 ^ farad 

P = .012 (P-P)/volt (rms) 

D . Charge per Droplet and Mobility 

From a total of 14 microscope photographs of droplets 
collected on a castor-oil-coated, 1/16 inch square, length 
of plexiglass, the most representative was selected for 
analysis and is shown in Figure 17 in Appendix C. About 6 
seconds were required to photograph the droplets after they 
were collected by being held for about one second, 10 inches 
from the 10-nozzle aerosol unit, at 40 psi air pressure. 

Due to the random nature of the droplet sizes, calculation 
of the charge distribution must involve rather gross assumptions 
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as to the most likely mechanism for actual droplet formation 
and charging. 

Droplets v/ere counted and measured over a representative 
section of Figure 17 comprising 1/8 of the photographed area. 
The results of the count are contained in table 5» , Appendix D. 

The calculations outlined in Appendix C were made for each 
of two different methods of distributing the measured output 
current with 300 volt (rms) input. From Figure 3* the peak 
current is about 1.85 pamps, with a flow rate of 1/3 cm^/sec 
corresponding to 40 psi air pressure from Figure 4b. 

Size-dependent Charge Distribution 

In order to make allowance for the likelihood of finding 
more charge on a large droplet than on a small one, average 
droplet diameter was calculated for three ranges of droplet 
sizes found in Figure 1 ?. Out of 60 "typical** droplets, 
the average diameter of the 6 largest was 22.5 pm. the next 
15 medium droplets 11.4 pm, and the average diameter of the 
39 droplets with diameters less than 10 pm was 5 pm. 

On the assumption that charge per drop is distributed in 
proportion to surface area the following charge per droplet 
and modified Rayleigh limit comparison results. 

(1) Large droplets (d-^22.5 pm) 

q-^~8.8 x 10 electrons/droplet 
q ,.5:3.35 x 10 electrons/droplet 
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(2) Medium dr oplets (d^js 11*4- pm) 

r 

q m ^ 2.28 :< 10 electrons/droplet 

CL 

q max » 1.21 x 10 electrons/droplet 

(3) Small droplets (d ?= 5 pm) 

s 

q 4.4 x 10^ electrons/droplet 
q max * J.5 x 10 5 electrons/ droplet 

Nu merical Average Droplet Diameter 

The usual convention for obtaining average droplet size 
in an aerosol is 



where 



N 

° 7 . d n /N 



(III.l) 



av S n=l 

IT = Total number of droplets in the sample 



Applying this definition to the data of Table 5 yields an 
average diameter 

d avg = 8.2.5 nm. 

the charge per droplet becomes 

q ~ 10^ electrons/droplet 
„ avg 



and 



q max = 7*5 x 10 ^ electrons/droplet 



Mobility 



Applying the foregoing results to equation (C.18) yields 
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an approximate average mobility, K (in units of coul-sec/kg 
m/sec 

or volt/m ^ * 

(1) L arge droplets (d^ = 22.5 pm) 

K = 5.62 x 10" 4 

(2) Mediu m dro plets (d = 11.4 pm) 

K = 1.85 x 10 “ 4 

(3) S mall droplets (d =5 pm) 

K = 8.1 x 10~ 5 

(4) Numerical average droplets (d ~ 8.25 pm) 

® * O 

K = 1.14 x 10~ 6 

E . Effe c t of Increasi n g Conductivity of the Water 

In view of the high rate of formation of the water droplets 
and a qualitatively observed fall-off of output current 
magnitude with increasing frequency, an observation of output 
current was made before and after a sudden increase in 
salinity of the water fed to the 10-nozzle aerosol unit. 

The apparatus was adjusted for an output current of 1.4 
pamp (P-P) with 40 psi air - pressure and siphon height of 
about 20 inches. 

A '‘normal" saline solution (.9 gm. NaCl per 100 cc. 

H 2 0) was substituted for the triple-distilled water (used 
in all other experiments of this report) and the output 
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current was observed to increase to 2.0 pamp (P-P). The 
current remained at this level as the system was being 
flushed with tap water. This experiment was repeated 
several times with consistent results. 



IV. DISCUSSION OF RESULTS 



Due to the extremely complicated nature of the actual 
interplay between system parameters, as previously discussed, 
the results must be viewed as more qualitative than quantitative 
The original objective was to investigate the synchronous ' 

EFD generator and, hopefully, to establish its feasibility as 
a practical energy conversion system. The results are 
encouraging in this regard, and are considered demonstrative 
of the fact such a system has sufficient potential to 

more than justify further development. 

The fact that the aerosol droplets are formed at a high 
rate from distilled water did not seriously limit the induced 
charge per droplet at the relatively low frequencies investi- 
gated. The simple instrumentation employed in measuring 
output current was not adequate for satisfactory determination 
of the response of the output current to increasing excitation 
frequency. This was due to the increasing capacitive 
admittance between the inducer electrodes and collector 
screen, and decreasing wave length of the charge stream as 
frequency increased. The former resulted in '’noise" pick-up 
swamping the charged aerosol signal, which was at the same 
time being reduced by the increasing width of the collector 
screen relative to the charge wave length. 

The 64 - 66 percent output/input current efficiency 
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achieved by this relatively simple generator design seems 
highly promising for eventual self-excited operation. It 
is quite likely that more precise construction and improved 
electrode design might increase this efficiency to 80 
percent or better. 

The efficacy of the induction method for charging the , 
aerosol droplets is considered to have been established 
by these results, at least up to frequencies of a few hundred 
cps. The field in the inducer section was 'not carried up to 
breakdown strength in order to avoid burning that capillary 
tube furthest off-center in its inducer. A machine-made 
nozzle array would possess a more uniform gap width and allow 
larger current output per nozzle. In addition, the output 
current can be increased by increasing the density of nozzle 
mountings, so that no further attempt was made to increase 
aerosol currents beyond those of Figure 3 for the purposes 
of this report. 

Another factor, equally as important as obtaining an 
increase in aerosol current, is the concurrent increase in 
droplet density within the channel. Since a practical energy 
converter will have to retrieve energy from the gas-dynamic 
flow, a high droplet density is required for effective 
electrical/mechanical coupling. The very low mobilities 
attained by these droplets, as described in the previous 
chapter, will only by useful if such is the case. 

In this study the aerosol was injected into a separate 
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air flow, but one obvious way to increase droplet density 
in the channel is to use only the air which comes from the 
charged spray apparatus. This might require changes in 
channel configuration for multistaging, but such modifications 
should be minor compared to the increased efficiency that 
would be realized. 

Another advantage of the sprayed aerosol method is the 
fact that recirculation in a closed system v/ould require 
no additional energy input. Since aerosol liquid is siphoned 
to the mixing nozzles, it need only be collected at the 
downstream end of the system and provided, v/ith a return path. 

In the previous chapter, two different methods were used 
in attempting to calculate the charge per aerosol droplet. 
While the size-dependent breakdown seemed intuitively to be 
a more accurate representation of the actual situation 
at the point of droplet formation, the resulting charge per 
droplet exceeded the modified Rayleigh limit used as an 
upper extremum. 

There are several possible explanations for this, in 
addition to the obvious conclusion that the assumption 
itself may be invalid. 

(1) The method used to collect and observe 
the droplets introduced a six second delay between the time 
of collection and that of photographing the results, 
the larger droplets it seems reasonable to assume that 



For 
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evaporation will just about off-set the flattening out of the 
droplets in the surface of the castor-oil-coating, and the 
observed diameter will be a fairly accurate measure of the 
original spherical diameter. For droplets on the order of 
5 microns and smaller, this will undoubtedly not be the case. 

In fact, it was observed that the lifetime of droplets 
on the order of one or two microns diameter, even on the 
castor-oil, was one or two seconds, at most. If a significant 
number of one to 3 or 4 - micron diameter droplets were lost 
in this way, it would explain the high calculated droplet 
charges, since the ratio of surface area to volume is inversely 
proportional to diameter. 

(2) Since the output current was observed to be 
somewhat conductivity limited, the actual circumstances may 
be controlled in part by the time it takes for a droplet 

to be formed. If all droplets are initially formed at the 
same rate, the result may be a tendency toward a more uniform 
value of charge per droplet, irrespective of size. 

(3) The fact that the charge per droplet 
calculated on the basis of an average size of 8.25 pm was 

an order of magnitude below Rayleigh's limit, may indicate that 
a truer mean value lies somewhere between the two assumptions, 

(4) From observation of the spray formation 
point it is apparent that the liquid emerging from a capillary 



I 
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tube is immediately drawn radially outward to meet the air 
jet. It may be that large droplets are formed before contact 
with the air jet, whereupon they are further atomized. In 
this event it would not be possible to achieve charge 
concentrations approaching Rayleigh's limit because of the 
increase in area to volume ratio through atomization. Here 
again, the actual mechanism may well be a random intermediate 
process where portions of the liquid ” sheet " are still 
intact upon reaching the air jet, and therefore some charge 
induction may take place within the actual atomization 
process . 

The interactions which take place in creating the aerosol 
spray are too complex to permit resolution of these questions 
here, but a more detailed investigation of the problem would 
be worthwhile if eventual optimization is carried out for 
this type of system. 

If there are a large number of very small droplets, 
as suggested in the foregoing, it would explain a good portion 
of the signal deterioration with distance down the conversion 
channel, described in Reference 19. As the droplets evaporate, 
the charges they carried become highly mobile, and very 
little loading ' would be required to pull them out to the 
channel walls. Since significant loading has not yet been 
investigated in the conversion channel, the experimental 
determination of mobility effects in the synchronous system 
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is recommended as an objective for further study. 

By increasing the operating frequency the efficiency 
of the system should be improved, since either the size 
of the generator can be decreased while retaining the same 
number of wave lengths, or more wave lengths can be accommo- 
dated and/or flow velocity increased. This would probably 
introduce no further constraints on the system other than 
a requirement for increasing the conductivity of the aerosol 
liquid. 

With size dependence a factor in the amount of charge 
carried per droplet, it can be seen, from the results in 
Chapter III, that for large aerosol droplet mobility becomes 
directly proportional to droplet radius. Optimization of 
the over-all system may involve compromise between increased 
current flow and maximum coupling v/ith the transport medium. 

In both cases, an increase in aerosol density is called for, 
and the optimum droplet size v/ill probably require a saturated 
transport medium to avoid evaporation in the conversion 
channel . 



V . CONCLUSIONS AND RE CO ’/.MEN DAT IONS 



A . Conclusions 

(1) The synchronous EFD energy conversion process is a 
feasible method of producing a-c electrical energy from 
fluid mechanical energy, 

(2) It may be more advisable to view the system as 
advantageous from the standpoints of reliability, light weight, 
and simplicity, rather than for expectation of high over-all 
energy conversion efficiency, 

(5) Achievement of electrical self-excitation should be 
feasible v/ith relatively minor improvement in efficiency 
of the experimental system, 

(4) The charge induction method used in this investiga- 
tion is capable of producing a distortionless sine v/ave 

of charged droplets without mechanical or electrical complexity 
or prohibitive input losses. 

(5) The use of distilled water did not impose a serious 
limitation on the magnitude of induced charge despite the 
high rate of formation of the aerosol droplets, 

(6) The charged droplet generator, as described, does 
seem to possess adequate potential for the further development 
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of a practical synchronous E?D system. 

( 7 ) The use of charged droplets larger than 50 to 100 
microns should probably be confined to ballistic systems. 

B. Recommendations 



(1) Significant loading of the charged droplet flow in 
the conversion channel should be investigated in order to 
better ascertain the nature of mobility effects and system 
performance capability for a synchronous EFD system. 

(2) System development seems to have reached the point 
where attainment of self-excitation is a feasible next step . 
This should most likely be attempted with increased nozzle 
density and elimination of the additional transport medium 
(flow fan) or the confining of its use to electrode boundary 
layer control. 

(3) Valuable information would result from detailed 
investigation of the electrical and atomization interactions 

at the nozzles.- The objective would be to determine the actual 
charge distribution over the range of droplet sizes found in 
the aerosol, to the eventual end of optimizing the design 
v/ith respect to mobility and current output. 

(4) Experiments . should be conducted with synchronous 
systems operating at higher frequencies and higher conductivity 
liquids. 
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(5) More sophisticated instrumentation for measuring 
aerosol current, especially at higher frequencies, should 
be developed. Refinement of the method of measuring droplet 
diameters should also accompany recommendation (3) above. 



VI. APPENDIX 



VI. APPENDIX 



A . Theoretical Background 

1 . Theory of the General S yn chronous Interaction 

To aid in analysis the model would be simplified to 
an effective thin stream of charge as in Figure 8. 




THEORETICAL MODEL for GENERAL 
SYNCHRONOUS INTERACTION 

The space between the charge wave and the channel walls 
is analagous to the air gap in a rotating electromagnetic 
synchronous machine. Just as the leading or lagging phase 
relationship between rotor and stator determines the direction 
and magnitude of energy transfer in the rotating machine, the 
same phase relationship exists between the traveling electric 
fields in the alternating EFD generator. This is intuitively 
apparent from Figure 8 and the mathematical analysis presented 
in the Appendix in Reference 19* A similar analysis of the 
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theory will oe outlined here for the general case as 
illustration of the interaction mechanism. 

The basic equations, in differential form, which must 
be satisfied within the channel are: 

Faraday's induction law V x E = 0 (A.l) 

Gauss* Law V • e o E = p = 0 (A. 2) 

These readily lead to the definition of a scalar potential 
which must satisfy LaPlaces equation between the channel 
walls. 

<P = -E (A. 3) 

whereupon : 

V 2 tp = 0 (A. 4) 

It is further assumed that the system can be considered 
semi-infinite in the Z-axis direction and that spacial 
derivatives with respect to Z are equal to zero. For 
simplicity, the symmetry of the problem will be used to 
justify solving the problem only in the upper half of the 
channel. The sinusoidal nature of the component waves will 
be observed by defining the charge and voltage waves as: 

a = Re[S 0 eJ (ut " 

j(o t - kx) -j 



and. 



<p = Re[$(y)e 



(A. 5) 
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Since the system is assumed, synchronous, the relationship 

u 

U_ = (A. 6) 

° k 

must be true. Solving LaPlace' s equation leads to 

^ A A 

<n(y) = A sinh ky + B cosh ky (A. 7) 

The boundary condition at the wall (y = d) requires that: 

V = A sinh kd + B cosh kd (A. 8) 

o 

The other boundary condition arises from applying Gauss 1 integral 
law to the upper half of the charge sheet at y = x = 0. 

That is: 



§ e Q E.da = J$p 0 dv = Jj — 2 da (A. 9) 

A 2 




— |1 O 

Figure 9. Gaussian Pillbox at charge sheet 



By symmetry with the other half of the channel, the only 
contribution to the closed surface integral is As o E(0) 
and, therefore: 
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A 

a 



E y (y=0) - 2 ~ 



6<p 

&y 



y=o 



(a. 10) 



It then follov/s that the other boundary condition requires 



A 




2s 0 k 



and, from (A. 8): 



B 



V o + 



a sinh kd 
o 

2 ke 0 



cosh kd 



(A. 11) 



(A. 12) 



Tn analyzing the forces on the charge shee.o, use will 
be made of the Maxwell stress tensor method as described in 
Appendix 1 in Reference 17. 

In this case the tensor at the charge sheet reduces to 
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(A. 13) 



SEE 
o x y 



I 2 < E y - E x> 



o 



\ 



0 

0 





+ 




where T represents force per unit area in the x-direction 
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on a surface normal to that direction, T. 



a stress in th 



xy 



x-direction on a surface normal to the y-axis, etc. 

The model for the differential section of half of the 
charge sheet 



shows that the forces cancel and the resultant stress 



traction in the v-direction, T . 

A* 

The symmetry of the model will require cancellation of 
all forces along the y-axis. However, in the actual system, 
where the width of the charge sheet is a significant fraction 
of the channel width, these forces result in detrimental 
spreading of the charges toward the walls. 

In accordance with (A. 3), 




o 



Figure 10 . STRESS TENSOR at CHARGE SHEET 




-E = Re 



jke^ [(-j)(A sinh ky + B cosh ky)(i x ) 



+ k(A cosh 




or: 




(A. 14) 
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Evaluating y = 0, 



E x = ♦ dkB^ ( “ 4 " kx > 

E y - - kAe^t ' kx) 

Nov;, the time average value of the traction becomes 

<t > = 1/2 Re[2e E e‘‘ ] 
x o x y 



(A. 15) 



(A. 16) 



Appropriate substitutions yield 



<t > = + ReCo 

x 2 



j8 »( 



A 



cosh kd 



. A A # 



0 Q tanh kd 
2ke_ 



■) 



(A. 17) 



Re Co o V ] 

2 cosh kd 



At this, point it becomes possible to readily observe the 
effect of the phase relationship between the voltage and 
charge waves. 

Rewriting (A. 17) in polar form: 

<T > k fs!s r c [>3(*A + 4>o - 

2 cosh kd 



where (•) and (•> are the phase angles of charge and 

C V 

voltage (atu t»kx=0), respectively. 



(A. 18) 



A further change of notation, together with multiplying 
<t x > by U , results in an expression for time-average power 
flow per unit area of charge sheet. 



.avg_ 



unit area 



ko V 

^ ooo f n e \ 

> = cos (— 5- - 6) 



2 cosh kd 



ka V 
0 0 0 

2 cosh kd 



sin6 



(A. 19) 



where 6 = (J) - ({) 



Taking the charge sheet as reference; 



6 = d> 

Y v 



It is now apparent that when the voltage leads the charge 
(6>0) by 90° , the force on the charge sheet is maximum in the 
direction of flow and the system acts as a pump. Similarly, 
when the voltage lags by 90 °, the energy transfer is greatest 
in retarding the flow and the system becomes a generator with 
flow dynamic energy of the charge sheet being converted to 
electrical energy in the walls. 

Figure 8 depicts the generating case. A more complete 
analysis of this case, the various modifications obtainable 
together with their further implications, and experimental 
results are to be found in Reference 19, a thesis in preparation 
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concurrently with this one. 

2 , Mobility Co n siderations 

The foregoing analysis is dependent on the assump- 
tion that the charges are rigidly coupled to the moving 
sheet. Evaluation of this assumption is extremely pertinent 
to a gas dynamic system since the coupling is not rigid, 
but is, infact, dependent on viscous slip forces which oppose 
the electrical retarding force. In this way the charged 
droplets interact with the gaseous medium to transfer energy 
either to or from the flow, as the case may be. 

It can nc., be ce'" r> that the amount of slippage required 
in order to produce the viscous force will be a loss in the 
net pov/er available from the system, since from equation (A. 19) 
power is proportional to the speed of flow of the charges. 

With the slip velocity defined as u, Mobility (K) has 
traditionally been defined as the ratio of slip velocity to 
electric field. 

K = — (A. 20) 

E 

Another detrimental effect of high mobility, is rapid 
lateral drift of charged droplets toward the side walls. 

For a droplet in the exact center of the channel, lateral 
forces would counterbalance, but obviously this would be an 
unstable condition even if it were physically, attainable. 
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Since mobility cannot be reduced to zero in a fluid system, 
lateral drift can, at best, be minimized or compensated, 
but probably never eliminated. 

Tbe mechanics of the interaction of ions and droplets 
in a fluid dynamic flow are well summarized in several of the 
references, particularly Reference 7 and Reference 15. 

The following is from a more detailed summary found in 
Reference 15. 

Since mobility is largely dependent on particle size, 
the usual three regions of categorization are by droplet 
diameter. 



(1) Small droplets (a <.002 microns) 

These include ions and other particles whose 
size is comparable to that of an air molecule. Here the 
coupling forces are predicted by molecular kinetic theory 
and both measured and theoretical mobilities decrease v/ith 
increasing pressure. The relationship for mobility derived 
from Longevin s theory is 



K « 3/4 



q 

6(8tiP^) 1//2 



(1 + m/Ll) 1/2 



where 6 = distance of closest approach 
m = ionic mass 
il = mass of gas molecule 
p = gas density 
P = gas pressure 



(A. 21) 
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V.'hile this equation predicts mobilities larger than those 
experimentally measured, the nobilities of ions are too large 
in gaseous media for adequate EFD efficiency. 



(2) Interme d iate droplets (.002 < a < .5 micron) 



Here a compromise between kinetic theory and 
viscous coupling is usually considered. The equation listed 
in both References 7 and 15, and derived from the Stoke- 
Cunningham law, is 

K = - (1 + — ) (A ??') 

6nv.a A 

where L = mean free path 

A = experimentally measured constant 
(on the order of .8) 

There. have been correction factors added to the second term 
for the larger values of L/A, but the equation as it stands 
is indicative of the tendency toward the Stokes law mobilities 
as droplet size increases. 



(3) Lar g e droplets (a>.5 pm) 

For large droplets Stokes lav/ leads to the simple 
expression for mobility, 



(A. 23) 



6 m p. a 
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For given charge per droplet it is evident that nobility can 
be decreased by increasing the size of the droplets. As 
previously mentioned, however, in the case of the induction 
system presented here the charge per droplet is also size 
dependent, as is the moisture accumulation on the channel 
walls. 

3 • Theoretical Limit on Gen e ral EFD Process with S li p 

A simple model of a direct EFD conversion section 
is shown in Figure 11. Irregardless of the actual interactions 
which may tab' 1 place within the section, use of Maxwell 
stress tensor analysis provides a convenient method of 
discerning the maximum energy conversion which could be 
realized. 

Several simplifying assumptions will be made, both for the 
sake of aiding analysis, and since they clearly provide a 
limiting case. 

(1) The electric field components in the x-axis 
direction are uniform across a constant area channel cross- 
section. 

(2) Energy changes in the fluid transport 
medium will be expressed as an equivalent pressure drop 
at constant flow velocity, U . 

(3) Flow frictional losses at walls and channel 
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ends will be neglected. 




o 



F igure 11 . GENERAL d-c EFD CONVERSION 
SECTION MODEL 

Y/hether the channel is rectangular or cylindrical, stress 
components .integrated over the side walls will cancel. The 
electrical stress magnitude on the channel section will be 
the difference between those on the surfaces normal to the 
x-axis. 



Equilibrium for the section requires the electrical stress 
and flow pressure differential to balance. 



2 

This expression will be maximum for E q = 0, and E^ at 
field breakdown strength, E^. 



T. 



x 




(A. 24) 



f_ (Ef - E^) = Pl - P 0 ^P 



(A. 25) 





2 



(A. 26) 
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The fluid dynamic power per unit area delivered to the 
conversion section will be 

4W u o * - 7 - 4 u o < A - 26 > 

while the electrical power density allowing for significant 
mobility of the charge carriei’s v/ill be 

L e * -J-4 (V 0 " V < A - 2? ) 

where L £ is the density of electrical power produced 
in the section and u^ is slip velocity at 
the inlet 

Substituting equation (A. 20) into (A. 27) 

' L e ' -7 E b (U o - K V ( A - 28 ) 

I 

where K v/ould be a measurable mobility according 
to the methods of the previous section* 

These results are obtained by Lawson v y ' in a way which 
demonstrates internal interaction of the system parameters, 
but the limiting case result is sufficient here to indicate 
the field strength limitation and effect of mobility in 
deterioration of efficiency. 
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B . Design of the a-c Charged Droplet Source 
1 • Theoretical Aspects 
Corona discharge 

The most satisfactory method, to date, for generating 
charged particles in direct current EFD generators, has been 
by condensation around, or direct convection of, ions formed 
from corona discharge at the entrance to the converter. The ' 
usual configuration provides an attractor electrode just 
downstream from a grounded needle (or needles). The field 
between attracror and n~edle ^ust be less th.qp breakdown 
but sufficiently high to .produce a space charge in the 
vicinity of the corona needle which is picked up by the 
v/orking medium and carried past the attractor. Since most, 
if not all, of the ions are not allowed to migrate to the 
attractor, the energy required to produce this current must 
be supplied by the medium. Figure 12 is from Reference 14- 
and shows the general arrangement of one type of corona 
configuration. 




Figure 12 . A CORONA DISCHARGE CONFIGURATION 
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C 14 ) 

Lawson' ' derives a theoretical relation for the 
magnitude of the current produced in a corona device which 
he then modifies to agree with experimental data. He begins 
by assuming a cylindrical space charge around the needle of 
radius R & , and constant charge density The additional 

assumption is made that the transport velocity is equal 
to half the flow velocity V. The current is then approximated, 



I = 



9 e 



V K R 



(B.l) 



and the field in the charge cloud (E q ) becomes 




(B.2) 



This, together with the definition of mobility, equation 
(A. 21), yields 



I 



£ # V Tt R, 



(B.3) 



Lawson obtained satisfactory agreement w ith experimentally 
obtained currents by including a factor of 30 percent. 

eV 2 x R c 

I = (.3) ~ — (B.4) 

K 



Substitution of suitable values corresponding to ions with 
V = 302 m/sec, R & = 1.59 mm, at an operating pressure of 
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19.4- atmospheres, yielded a value of I' equal to 120 pa. 

In turning his attention to colloidal charge carriers, 
Lawson assumed slip to he negligible (V = V) and E was 
therefore allowed to approach breakdown, E^. Using the same 
conditions as before, and dividing the resulting expression 
for I c0 n. 0 ^£ By that for 1^ , Lawson obtains the ratio, 



T col 

^ion 



= 1080 /V 



where V is expressed in m/sec 



(B.5) 



His experimental values were found to be in general 
agreement and colloid currents more than adequate for a 
generator contemplated under the limitations imposed by space 
charge effects within the channel^“’ ^ « 

Some kind of modified corona method may eventually 
prove feasible for alternating polarity ion generation, but 
at the state of present technology, the much lower breakdown 
potential for positive ions seem3 to be a prohibitive 
limitation^ 



Induction method 

It was decided, for this investigation, to use capacitive 
coupling betv/een an inducing electrode and sprayed water, 
droplets to induce charge of opposite polarity on the droplets 
as they break av/ay . A considerable amount of research has been 
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conducted in the study of electrical spraying of conducting 
liquids (References 1, 5, 9* 10, 16) but this has also been 
with d-c fields and largely for application in thruster 
design for space vehicles. Droplets formed due to electrical 
stress instability could be expected to produce a severely 
non-uniform charge distribution if the signal on the inducing 
electrode were sinusoidal rather than constant potential. 
Since the aim of this research v/as to produce a harmonic 
free, " clean*' sinusoidal charge wave, the instabilities were 
avoided. 

The limiting value of charge that a liquid droplet can 
support is a well known expression derived by Rayleigh^^. 

It is the limiting surface charge density just short of 
overcoming the surface tension holding the droplet in 
equilibrium. 

The expression for Rayleigh's limiting value of charge 
on a spherical droplet is 



q = 8 tc e 0 y (B.6) 

A similar limiting expression can be readily obtained for the 
case of a droplet being formed under the combined stresses 
of electrical field and internal pressure to overcome surface 
tension. A simple model is shown in Figure 13 which uses 
the rather gross assumptions of hemispherical shape as the 
droplet is released, and a uniform electrical field acting 
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normal to the liquid surface. 



Y p Y 




Figure 15 . MODEL for MODIFIED 
RAYLEIGH LIMIT 

E 2 

The electrical stress on the surface is — E^ and the 

2 

equation of o aui librium for the hemisphere just prior to 
separation is 

*/2 e 

y2Tca' = 7t(a' )^P + J — ~ E^cos© 2ma sin©(a' d©) 

0 2 

,2 



= n(a' )^[P + 



e E 



(B.7) 



If the liquid is considered to be a "perfect *' conductor, 

Gauss' lav; at the surface of the drop requires that electrical 
displacement (e o E) be equal to the surface charge density (a). 
Equation (B.7) becomes, 



2y = a' [P + 



2e. 



(B.8) 



(q/2ma’ J ) 2 

= a' [P + ] 

2e 
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After separation the dropl .t radius , a, is obtained by equating 
the hemispherical volume to that of the spheroid, whence 



Substituting (B.9) into (B.8) and solving for the charge 
per droplet 



It can be seen that for this model, if the pressure differen- 
tial is neglected, the maximum charge per aropret becomes 
approximately 70 percent of Rayleigh's limiting value. 7/hile 
this model could not be expected to accurately represent the 
complex situation of the aerosol spray, equation (B.10) will 
be used as a basis for comparison since it would undoubtedly 
provide an upper limit, 

2. Aerosol G e nerator Design 

The two nozzle designs are similar enough that only 
the 20-nozzle unit will be described in entirety. The 
electrode configuration of the ''fine droplet" unit is shown 
in Figure 16. 

Figure 14 is a schematic of the over-all configuration 




(B.9) 



q = 4nf2?fe 0 



( r 



2 2/3 



(B.10) 



of the charged aerosol generator. The major elements are 
the compressed air harness, the v/ater siphon line and nozzle 
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supply headers, and the two 5-nozzle spray units. 

The channel section formed by the generator unit is 5 
inches long and tapered in order to adapt the flow fan 
section to the main conversion channel, and to provide some 
convergence of the air flow as it picks up the aerosol. 

The entrance cross-sectional dimensions are 6 inches by 
2— 5/4 inches and the exit dimensions are 6 inches by 2-1/4 
inches to match those of the downstream conversion section. 

The general simplicity of the design is perhaps best 
illustrated by the fact that all parts except tubing were 
constructed of plexiglass, with a drill press being the 
only power tool required. 

Figure 15 shows the detailed construction parameters 
of the 20-nozzle aerosol unit. As previously observed, the 
peculiarities of the over-all system provided the majority 
of the design constraints. The only basic requirements of 
the aerosol generator itself are the positioning of a water 
carrying capillary tube in the center of the cylindrical, 
high velocity air jet, and the arrangement of the charge 
inducing electrode in close proximity to the mixing point. 

To produce a uniformly fine spray, the water must be 
siphoned to the nozzle by the concentric air jet. In general, 
the greater the siphon distance the finer the spray. The 
siphon height for ths charged aerosol generator also has a 
small effect on the current output, since the flor/ rate is 




water 



tEfc AIR 



F-'CV-'o 14. SCHEMATIC of 
OVER-ALL AEROSOL GENERATOR 



AIR 





Fi gure 16 . 

CONFIGURATION of 10-NOZZLE 
INDUCER 



Fi gure 15 



SCHEMATIC of 20 -NOZZLE UNIT 
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greatest for small siphon distances, A siphon height of 
about 23 inches seemed to be the best compromise for this 
system. 

The size of the air feed channel in the nozzle assembly 
has to be large enough to ensure a reasonably uniform air 
pressure feeding all of the individual nozzles. This was 
especially true with the limited laboratory compressed air 
supply. 

3 • Measuring; Drop l et Size 

In order to calculate the average cnarge per droplet and 
charge-to-mass ratio for the traveling wave particles, it is 
necessary to know the average size of the aerosol droplets. 
Probably the simplest way to accomplish this, in theory, is 
to collect the droplets on a slide and photograph them under a 
microscope. This is the method that was used in obtaining the 
data for this thesis report. 

When highly accurate results are required, the process 
of collecting the droplets becomes critical. Under such 
conditions, the aerosol is usually allowed to fill a saturated 
container and a very srender collecting glass is passed 
carefully through the fog at an accurately controlled rate. 

If the collector presents an excessively large cross-section 
in the direction of notion the smaller droplets are swept 
aside and only the larger ones penetrate the boundary layer 



to strike the slide 
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Another problem which must be faced, irregardless of the 
accuracy required, is the short lifetime of micron sized 
droplets in "dry" air. In order to prevent rapid evaporation 
the slide is coated with castor oil which has been in contact 
with water long enough to have become saturated. 

In view of the fact that the aerosol is injected into a 
dry air stream when used in the present synchronous system, 
the droplets were collected in the dry air environment on 
the coated slide and the results are considered more or less 
representative of the conditions as the aerosol enters the 
converter section. Since the droplets require only 1/20 sec, 
to traverse the channel, it was assumed that further 
evaporation within the conversion section could be neglected. 
Figure 17 shows a typical microscopic view of the collected 
droplets at about 330 X* 

C. Calculations 

1 . Electrical Capacitance of Aerosol Generator 

The electrical circuit model of the apparatus used to 

measure the equivalent capacitance is shown in Figure 18, 

where V was a variable frequency (10 - 1200 cps) sinusoidal 
s 

voltage source. V.^ and V 2 were fed to a dual trace oscillo- 



scope 
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Figure 16 . CIRCUIT MODEL of INPUT TEST APPARATUS 



Frequency Response 

In order to construct a Bode diagram with its break- 
point within the available frequency range , R was set a 48 
megohms. In LaPlace transform notation, the admittance 
relation I/V g is 



I 2. 10' Cs 
V 23 V, RCs + 1 

S 1 



(c.i) 



In terms of magnitude alone for s = joj, 

- 1 1 1 = |Vsl CH 2 oV ; uV2 



(c.2) 



taking the logarithm of both sides and multiplying by 10, 



(C.3) 



10 log 1Q }l| = Constant + 10 log^to .. 51og 1( - ) [(RCto) + 1] 
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The constant can be ignored in analyzing the data, 
while the remaining tv/o terms should provide tv/o asymptotic 
regions. At very low frequencies a plot of 10 log-^Qjll vs. 
10 1 °S;lo w should approach a straight line of slope 10 
decilogs per decade of frequency. As RC OJ approaches 1 the 
curve should begin to " break ", and at higher frequencies 
approach a constant value. This can be seen from the limit 
of equation (C.2) as to goes to infinity. 

The phase plot of current lead angle vs. 10 log-^Qt-o 
provides an additional check on the preceding, because the 
break point at RCto= l also means a 4-5 degree phase angle 
for the denominator term of equation (C.l). The results are 
plotted in Figure 6 and the data summary is contained in 
Appendix D. 

Numerical Check 

Equation (C.2) was used to directly compute the value of 
capacitance (C) for several values of resistance (R) at a 
constant frequency (f) of 125 cps. 

Letting R l = (R + .5) meg- ohms, and rewriting equation 
(C.2), 

2 V 2 C2mf 

25 V 1 .10 6 [(R‘C27tf) 2 + 1] 1/2 



(C.4) 
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then 




(c.5) 



Substituting f = 125 cps and holding constant at 

10 volts 



C 



(i2.7)v 2 .io~ 12 




(C.6) 



I 

With R = 10.2 x 10 , the measured value of V 2 was .62 volts. 
Then, 



7.88. 10" 12 

G s: 

1.15 



or 



C » 6.85* 10"- 1 - farad 
Theoretical Model Capacitance 

As an order of magnitude check on the measured capacitances, 
a crude model was calculated. The capacitance between an 
individual nozzle and its electrode was compared to a 
cylindrical capacitance of similar dimensions neglecting 
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fringing, as shown in Figure 19 below. 



h — i — -H 




Figure 19 . CYLINDRICAL CAPACITANCE MODEL 



From the fact that the field in cylindrical coordinates is 
inversely proportional to radius and from application of 
boundary conditions at r = R^ and R , it readily follows that, 

e 2x1 

c = _o (C.7) 

lnCR^) 



For both nozzle assemblies the ratio R q /R^ is approximately 
2. Substituting numerical values into (C.7) 



C ' 



(1/36tc) •10~9(f/m)«27il(m) 
In 2 



(,&) 10 _1 ^l(cm) farad (C.8) 

where C* is the capacitance of an individual nozzle 

For the 10-nozzle unit the order of magnitude of capacitance 
will be 

C cs 8* 10”^l(cm) farad 



(C.9) 
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2 . Calculation of Charge -per Dr oplet 




As a further refinement of the base of comparison 
represented by equation (B.10), a minimum value of pressure (P) 
(after allowing for a siphon height of 2 feet) will be 
presumed. For an air inlet pressure to the 10-unit nozzle 
of 40 psi., siphon height in excess of 23 inches was observed 

•7 2 

to be at least 1/14 atmosphere, or 73*05 x 10^ dyne/ cm - 
pressure. For simplicity of calculation in equation (B.10) 
therefore, this wi 1 l be used as reference pressure. 



Substituting the appropriate constants into equation 
(B.10) 





where "a"is droplet radius in centimeters 



Expressing’^" in microns (10" /+ cm), the expression becomes: 
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(C.ll) 



Equation (C,ll) represents the extreme upper limit from 
Appendix B (P= 0) by neglecting the negative term under the 
square root , or could be adapted to other pressures by 
multiplying that term by 13*9 P, where P is expressed in 
standard atmospheres. 

Droplet Flow Rat e 

The measured flow rates were in terms of cm^/sec which 
can be converted to gallons/hour through the conversion 
factor 



count droplets by size and use a straight numerical average 
as the mean droplet diameter. On that basis the data taken 
from Figure 1? results in an average value of 



Q(GPH) = Q(cm 5 /sec).(.95) 



(C12) 



Calcula t ion of Charge per Dropl et 



In analyzing "fog* 1 droplet data, usual practice is to 



135 + 171 + 189 



= 8.25 pm 



60 



For the purposes of calculating charge per droplet, 



this is not necessarily the most meaningful type of average. 
In order to investigate a size dependent measure of the 
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charge per droplet, the distribution will be calculated on 
the basis of average large droplet”, 'medium droplet", and 
average "small droplet". 

The average large droplet diameter is s & = ^^*22.5 pm, 

the average medium diameter a d^ = * 11.4 pm, and the 

average small droplet diameter = & = --§§ 5 pm. The total 

s 2 ? 

number of droplets per sec. for a given flow rate is 
defined, 



N S (6 + 15 + 39)x = 



Q 



6 CG ^1 



I5d^ + 39 O 

m s 



(C.13) 



Using the flow rate for the 10-nozzle unit with air 
pressure of 40 psi,from Figure 4. 

12 

Q = . v cmVsec = — (pm)^/sec 
3 3 

Equation (C.13) yields 



6-10 



12 



x = 



60(3)tc[6(22,5) 5 + 15(11. 4) 5 + 39(5) 3 ] 



= 1 . 11 * 10 ^ (60-droplet units/sec) 



Since the current per droplet will be roughly proportional 
to surface area, the output current must be apportioned 
between the 3 droplet: classes according to their respective 
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fraction of total surface area. 



Total surface area = 



A - nx[6d^ + 15d^ + 39dg] 



(C.14) 



For the case in point, 



A = 5965 -nx (pci ) 



Large droplet fraction 5 f 



f. = - 



6dJ 

5965 



(C.l^ 



= .509 



The medium droplet fraction is 



f m * 



1950 

5965 



= .327 



and the small droplet fraction becomes 



f s = .164 

From Figure 3* the output current for a 300 volt input v/ith 
the 10-nozzle unit, is I = 1.85 U a i peak value. 



Large Droplet Charge 

The fraction of output current carried by large droplets, 
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>6 coul 



2 I ol = 1.85-10“° (.509) =(.94)10 



-6 coul 
sec 



The charge per droplet becomes 



q l = 



Id coul/sec 
6 x droplets/scc 



(C.16) 



= 1.41« 10 coul/droplet 



Since the electronic unit of charge is e = 1.602. lO^^coul , 
the approximate number of electrons per "large” droplet is 



ul_ ~ 2.2 !0 6 ( ‘-tectrons/drop'i 



Medium Droolet Charge 



I = .606*10"" coul/sec 
on 



and equation (0.16) becomes 






— 1 - — = 5.65-10”^^ coul/droplet 



15 x 



and 



g 2.28*10^ electrons/droplet 



Small Droolet Charge 
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I os = .30r*10~° coul/sec 



Q = — ^ — = 7 . 03 - 10” 14 coul/droplet 
s 



39 x 



and 



< 1 - ^4.4.10^ electron/droplet 



** Modified" Rayleigh Limit 

Using eolation (C.ll), and neglecting the pressure term 
(1) Large droplets 



q max js 14.2. icr 15 (n.25 ) J/2 

= 5.37*10 -1 ^ coul/droplet 
= 3.35‘10 electrons/droplet 



(2) Medium droolcts 

q max « 14.2.10- 15 (5.7) 3/2 

c 1.9^-. 10~ 15 coul/droplet 
= 1.21*10^ electrons/droplet 
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(3) Small droplets 

W ~ W-2'10" 15 (2.5) 5/2 

= 5. 62 •lO”^’ coul/droplet 
= 3»5*10^ electrons/droplet 

(4) Numerical " Average " Droplets 

<W = W-2-10- 15 (4.1J) J/2 

~ 1.2»iO"''*' 3 coul/droplet 
= 7. 5* 10 3 electrons/droplet 



Charge per Drorlct for Numerical 'Average" Droplets 



Droplet flow rate can be expressed: 



Q,(droplets/sec) = 



10 12 Q 

4tc a 3 

3 



CC.17) 



_ 10 12 (l/3) 

•;4n (4. 13 ) 3 
3 

= 1 . 13 * 10 ^ droplets/scc 



charge per drop becomes > 
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1.85*10“^ coul/sec 

- 

1 . 15*10 droplets/sec 
~ 1.65 * 10*”^^ coul/droplet 
~ 10 4 electron/droplet 



Mobility 

Substituting numerical values into equation (A, 25) 



K = 



2q (coul) 
d(pm) 



(■i) 



184*10 



-6 



Km 



sec cm 



10 “^cm 

pm 



10“ 5 Kg 

gm 



= Q ( coul) 5.77-10 9 
d Cum) 



(C.18) 



D. Summary of Data 



TABLE 1. 



(Rms) 



Single Commercial Nozzle 

Air Flow 

Press Siphon Height Rate 

(psi) (GPH) 



(ua, P-P) 



100 


47 


24 


.27 


.13 


150 


52 


24 


.28 


.19 


200 


51 


24 


.28 


.24 


250 


50 


24 


.28 


.3 


500 


50 


24 


.28 


.37 


200 


75 


24 


• 56 


.28 


200 


70 


24 


.54 


.26 


200 


65 


24 


.32 


.25 


200 


60 


24 


.30 


.25 



V 1 

80 

100 

150 

200 

250 

300 

100 

100 

100 

100 

100 

100 

100 

100 

ICO 

100 

I 
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TABLE 2 . 



10-Unit Nozzle 



Air # 

Press 2nd I 

0 



Flow Rate j 

(cn , /sec) x o 



40 


.75 






40 


1.25 




.8 


40 


1.75 


.333 




40 


2.5 




1.4 


40 


3.0 






40 


3.8 






50 


1.3 


.358 


.83 


45 


1.3 


.37 


.83 


40 


1.25 


.333 


.81 


35 


1.15 


.233 




30 


.9 


.20 


.69 


25 


.85 


.143 




20 


.7 




.5 


15 


.44 




.23 


55 






.81 


60 


■vvv 

( electrode 


radius reduced 


.79 



by about 1/3) 



TABLE 3 . 



Impedance Frequency response for 10-nozzle unit 
(V^ = 9.4 (P -P) Refer Figure 18) 



Freq ^2 Phase Angle 

(cps) (P -P) (V 2 leading V x ) 



40 


.2 


70 


80 


.36 


63 


100 


.43 


60 


150 


.60 


60 


200 


.74 


57 


300 


1.00 


54 


400 


1.20 


46 


500 


1.36 


42 


600 


1.44 


38 


700 


1.52 


35 


800 


1.6 


32 


900 


1.68 


28 


1000 


1.8 


23 


1100 


1.86 


20 


1200 


1.9 


18 



